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Abstract 

The first run of the LHC was successful in that it saw the discovery of the elusive Higgs boson, a particle that is consistent 
with the SM hypothesis. There are a number of excesses in Run 1 ATLAS and CMS results which can be interpreted 
as being due to the existence of another heavier scalar particle. This particle has decay modes which we have studied 
using LHC Run 1 data. Using a minimalistic model, we can predict the kinematics of these final states and compare 
the prediction against data directly. A statistical combination of these results shows that a best fit point is found for 
a heavy scalar having a mass of 272^g 2 GeV. This result has been quantified as a three sigma effect, based on analyses 
which are not necessarily optimized for the search of a heavy scalar. The smoking guns for the discovery of this new 
heavy scalar and the prospects for Run 2 are discussed. 
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1. Introduction 

With the discovery of a new scalar boson (which will be 
denoted by h ) at the Large Hadron Collider (LHC) [1, 2], 
new tasks and explorations have come to the fore. Both 
the ATLAS and CMS experiments have been keenly ex¬ 
ploring the properties of the new scalar, and will continue 
to do so in the high-energy run(s). This includes tests 
on various coupling strengths as well as spin-CP proper¬ 
ties. Overall, the global fits of the data indicate that the 
properties are mostly consistent with what the Standard 
Model (SM) of particle physics expects for a Higgs boson. 
Nonetheless, there are certain features and excesses in the 
available data that warrant detailed attention. Table 1 
summarizes the measurements made by the ATLAS and 
CMS experiments that are considered here. These results 
comprise the measurements of the differential Higgs bo¬ 
son transverse momentum {pth) (where h is the SM-like 
Higgs.), searches for a di-Higgs boson resonance, the Higgs 
boson in association with top quarks, and VV resonances 
(where V = Z, W ± ). These final states are considered 
against the hypothesis of a heavy scalar boson {H). 
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In order to describe the shape of the p^u distribution, 
it is necessary to introduce decays in which at least one h is 
produced: H -A hh, h\X where \ would be a dark matter 
candidate, leading to the production of the Higgs boson 
in association with missing energy. The latter could be 
realized through the decay of some intermediate particle. 
This hypothetical intermediate particle (the existence or 
nature of which we make no statement about) could also 
decay into a pair of light hadronic jets. In either case, a 
distortion of the pxh spectrum would be expected. 

Both ATLAS and CMS have observed excesses in the 
production of a SM Higgs boson in association with top 
quarks. The new boson would naturally be produced in 
association with top quarks. In addition, with a small 
branching ratio to VV, the effect of negative interference 
in single top associated production is suppressed, and we 
would expect a large cross section for a heavy scalar being 
produced in association with top quarks. These effects 
would yield an explanation for the excesses which are seen 
in the data. 

While none of the excesses summarized in Table 1 is 
significant enough on its own, it is tantalizing to study 
their compatibility as a whole with the decay of one new 
heavy boson. In fact, viewed under this hypothesis, the 
overall size of the excess could be significant. We adopt 
a bottom-up approach and wait for an overseeing theory 
until the observed results are further consolidated. The 
fit of the available data presented here by us will hope¬ 
fully enable people to put all pieces of the jigsaw puzzle 
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Result 

Publication 

Differential Higgs 
boson pt spectra 

ATLAS 

Fiducial cross section measurements on h —► 77 [3] and h —¥ ZZ* — > [4] 

CMS 

Fiducial cross section measurements on h —> 77 [5] and h —► ZZ* —» [ 6 ] 

Di-Higgs boson 
resonance searches 

ATLAS 

Limits on H —> hh —> bbTT, 77 WW*, 77 bb and bbbb [7] 

CMS 

Limits on H —>■ hh —)• 77 bb [ 8 ], bbrr [9] and multi-lepton [10] 

Top associated Higgs 
boson production 

ATLAS 

Limits on h —> 77 [11] 

Measurements on multi-lepton decay channels [12] and h —»• bb [13] 

CMS 

Measurements on h —> 77 , h —> bb and multi-lepton decay channels [14] 

H —> VV decays 

ATLAS 

Limits on H —> WW [15] and ZZ [16] 

CMS 

Limits on H —> WW and ZZ [17] 


Table 1: A list of the experimental results which were used to help constrain the relevant parameters of the proposed model. In the interest 
of being as unbiased as possible, these results were selected regardless of whether they hint at physics beyond the SM. 


together once a larger volume of results accumulates. It 
is important to note that some of the experimental anal¬ 
yses used in this study were not optimized for the search 
of a heavy scalar boson. Therefore, the sensitivity to the 
search is not maximized, rendering conservative the results 
reported here. 

The salient features of the suggested scenario are sum¬ 
marised in Section 2. In Sections 3 and 4 the main tools of 
our analysis are discussed with Section 3 focussed on the 
statistical machinery used. The results are presented and 
discussed in Section 5 including a discussion on smoking 
guns and prospects for Run 2, followed by a brief conclu¬ 
sion in Section 6 . 


2. A suggested scenario 

The features and excesses which are seen in the data 
have been treated in such a way that they could be ex¬ 
plained purely by physics beyond the SM (BSM). We pro¬ 
pose a scenario which is a BSM extension to the SM, al¬ 
lowing us to write 


-S? = -£?sm + -2bsm, (1) 

where all of the new interactions and states are encoded 
in 2 z?bsm- 

The simplest approach is to treat the new interactions 
as arising from effective couplings. We have assumed h 
to have SM interactions with fermions and gauge bosons. 
The sectors of the proposed BSM Lagrangian involving the 
new scalars (omitting the usual mass and kinetic energy 
terms) include, 

-^bsm 7) + -%t + -^Qj (2) 

where terms _2y, -2 t an d 2 S?q are the Yukawa, trilinear and 
quartic interactions relevant for this study, respectively. 
These sectors are defined as follows, 

*h = -\ P g *Z G^H + P v kZ v V^H, (3) 
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hhxx - yA a „ HH X X, 


( 4 ) 

( 5 ) 

( 6 ) 


where H and % denote the heavy scalar and the DM can¬ 
didate respectively (the latter is assumed to be a scalar for 
illustration), and v = 246 GeV is the vacuum expectation 
value that is responsible for the W- and Z-boson masses. 
This can be looked upon as a variant of Higgs boson portal 
scenarios [18, 19, 20]. This Lagrangian could in principle 
emerge as an effective theory after electroweak symme¬ 
try breaking in any gauge-invariant extended scalar sector. 
The second term in Eq. (3) is summed over the weak vec¬ 
tor bosons Z and W , and the k factors are the SM-like 
couplings, with = a s /(3i tv) and kZv — 

We deliberately make no statement about the gauge 
quantum numbers carried by H, but just postulate that 
the above terms stay after electroweak symmetry break¬ 
ing. We set Pgi^Z, be the strength of the effective 
gluon-gluon coupling of H. In situations where there are 
no additional effects over and above the top-mediated tri¬ 
angle diagrams contributing to this effective interaction, 
Pg = VitH/Utth w here y tth is the SM top Yukawa cou¬ 
pling. There would also be a similar relation for the bot¬ 
tom Yukawa coupling y bbH , but this has been counted as 
negligible since the effect of bottom quarks in gluon fu¬ 
sion loops is small. The production of H is made to occur 
through gluon fusion and its rate can therefore be con¬ 
trolled by varying /3 g . Likewise, the HVV couplings can 
be tuned by varying /3 V . 

It should again be remembered that we are not making 
any definite statement on the origin of each term. Thus, 
as we shall see below, the Hhxx (effective) coupling is re¬ 
quired to be on the high side, keeping the observed data 
in mind. We do not rule out the possibility of this be¬ 
ing due to the participation of some real particle in the 
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intermediate state, as discussed above. 

Major constraints on the Lagrangian parameters stem 
from the observations of the relic density of DM [21] and 
the DM-nuclei inelastic scattering cross sections [22], These 
constraints are controlled by the two model parameters 
m x and A h xx - It is found that both of these DM con¬ 
straints can simultaneously be satisfied for a narrow choice 
of the parameters m x ~ [55 — 60] GeV for very small 
A h X x ~ [0.0006 — 0.006] [23]. This keeps the invisible decay 
width of h well within the observed limits. Other model 
couplings remain unconstrained by these observations. 

Within this simplistic framework, one expects the pro¬ 
cess pp —>• H —>■ /iyy t° generate an enhanced pr, owing to 
the fact that h now recoils against a pair of invisible parti¬ 
cles. The presence of a H\X coupling opens the potential 
for detecting invisible decays following the methodology 
suggested in Ref. [24] . 

3. Statistical formalism 

The BSM prediction constructed using the proposed 
model was fit against four classes of constraints. The con¬ 
straints considered are differential Higgs boson px spectra, 
the limits on di-Higgs boson production through a reso¬ 
nance, various limits and measurements on top associated 
Higgs boson production, and the limits on a heavy scalar 
decaying to vector bosons. For each class of constraints, 
results from both ATLAS and CMS were used to avoid 
bias. These results are summarized in Table 1. A simulta¬ 
neous fit was done in terms of calculating and minimising 
a combined x 2 value while varying /3 g for different to# 
hypotheses. Statistically, two types of results were dealt 
with: measurements and limits. 

A calculation of y 2 for a measurement is straightfor¬ 
ward. Given a measurement p with its associated error 
A/i, one can construct a y 2 by testing the measurement 
against a theoretical prediction and its associated theoret¬ 
ical uncertainty, given as p th and Ap th , respectively. The 
uncertainties from the measurement and the theoretical 
prediction are assumed to be independent and are added 
up in quadrature, allowing us to calculate the y 2 as 

2 = (/i — /x th ) 2 

* (Ap) 2 + {Ap th ) 2 ' {> 

To calculate a y 2 from a result in the form of a 95% con¬ 
fidence limit (CL), we need only assume that given some 
measurement p with its expected and observed limits p exp 
and p ohs respectively, the y 2 is Gaussian in p. Then, as¬ 
suming the null hypothesis has p = 0, we can extract the 
mean of the distribution as p ohs — p exp . This is treated as 
the excess in p which can be tested against a theoretical 
prediction p th . Its error is inserted as A p = p exp / 1.96, 
where the 1.96 arises from the fact that 95% confidence 
corresponds to 1.96 units of standard deviation in a Gaus¬ 
sian distribution - this approach is used in Ref. [25] and 


other references therein. Using this, the y 2 is calculated 
as 


X 


2 


(p° hs - p exp - p th ) 2 

(p exp /lM) 2 


( 8 ) 


Using the definitions in Eqs. (7) and ( 8 ), a combined 
y 2 was constructed by adding up the contributions from 
all of the results presented in Table 1. This procedure is 
described in the following section. 


4. Methodology and tools 

The minimal model we have described was built first 
using FeynRules [26] and then passed to the Universal 
FeynRules Output [27] such that event generation could 
be performed at leading order (LO) in MadGraph5 [28]. 
Computations relating to DM constraints in the model 
were carried out using mi crOMEGAs [29]. 

Fitting the ATLAS and CMS Higgs boson pt spec¬ 
tra with the BSM prediction was accomplished as follows. 
Events with the ft.yy and hh final states were generated 
from pp collisions through an H s-channel in MadGraph5 
and showered appropriately using PYTHIA 8.2 [30]. Since 
the ATLAS [3, 4] and CMS [5, 6 ] Higgs boson pt spec¬ 
tra were constructed from fiducial volumes of phase space, 
it was important that BSM prediction went through the 
same event selection. This was done using the Rivet [31] 
analysis framework. The total LO cross section of the BSM 
prediction was enhanced to NNLL+NLO (next-to-next-to 
leading log plus next-to leading order accuracy) through 
multiplication by an appropriately calculated k-factor de¬ 
termined from Ref. [32], 

The BSM prediction was tested against the SM-only 
prediction. For the SM, fiducial acceptance factors for 
gluon fusion (ggF) are given in Refs. [3, 4, 5, 6 ]. The ggF 
Higgs boson px prediction was generated at NLO using 
MiNLO HJ [33]. This prediction was further reweighted 
to NNLO accuracy and scaled to the NNLL+NLO cross 
section from Ref. [32]. The other less prominent Higgs 
boson production modes (VBF, Vh and tth , which are 
collectively referred to as Xh ) were taken directly from 
the ATLAS and CMS publications. 

The bulk of the production cross section is in the in¬ 
termediate pxh range where ggF is the dominant produc¬ 
tion mode. In this mechanism, QCD radiative corrections 
play a critical role in generating pxh- The Monte Carlo 
(MC) used to simulate ggF describes the pxh distribution 
at NNLL+NLO. Recent results on the NNLO corrections 
on ggF+lj production indicate that, although moderate, 
corrections are still significant [34, 35]. NNLO corrections 
with respect to NLO can be as large at 25% in the range 
of interest. In order to accommodate these corrections, a 
conservative approach is implemented. The pxh distribu¬ 
tion with pxh > 30 GeV is corrected with the NNLO/NLO 
k-factors provided in Ref. [34]. The MC described is nor¬ 
malized to the total ggF cross section at NNLO. The first 
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Figure 1: Fits to the fiducial differential distributions of the Higgs boson transverse momentum using the ATLAS diphoton [3] (left), the 
ATLAS h —>■ ZZ* —> [4] (top right), the CMS diphoton [5] (middle right) and CMS h —> ZZ* —> 4 i [6] (bottom right) decays (see text for 

detailed description). The mass points considered here are the best fit values of m\ = 60 GeV and m# = 270 GeV. In the interest of saving 
space, only one full plot has been shown, with the other three are plots of a ratio between the data and BSM prediction. 


complete calculation of the total ggF cross section at N 3 LO 
is now available and it is indicative of small N 3 LO/NNLO 
k-factors and scale variations [36] . For this reason the cross 
section with prh < 30 GeV is re-scaled appropriately so 
that the total cross section does not exceed the total cross 
section by more than 2% with respect to the calculation 
at NNLO. The scale uncertainties assumed in this analysis 
remain at NNLO for the total cross section and at NLO 
for the pTh, while the PDF uncertainties were conserva¬ 
tively taken from Ref. [32] . Other production mechanisms 
of the SM Higgs boson do not play a critical role in the 
region of the phase-space under study, and were therefore 
left unmodified. A similar approach is taken for the BSM 
prediction. The pr of the heavy scalar is reweighted by 
determining a reweighting function from an NLO calcula¬ 
tion using OneLOoP [37] in MG5_aMC@NL0 [38] compared to 
the result obtained with the shower. Overall, the shower 
does a reasonable job, matching the LO prediction within 
20%. The effect of the these corrections on the transverse 
momentum of the Higgs boson from the decay of H is 
small and it reduces to a positive shift of about 3 GeV. It 
is, however, important to note that the jet multiplicity of 
the H boson in this setup is significantly larger than that 
characteristic to h. This implies a significant reduction of 
the jet veto survival probability. 

Excesses in top associated Higgs boson production were 
also included in the fits. Associated th production in the 


SM is suppressed due to the negative interference induced 
by the relative sign of the Yukawa and hWW couplings 
(see Ref. [39] and other references therein). If the HWW 
coupling is relatively suppressed, this negative interference 
is reduced, so that its cross section becomes comparable 
to that of ttH production. For this reason, /3 V was set 
to a small value (order of 10 —3 ) and tH cross sections 
were determined at LO in MadGraph5. These cross sec¬ 
tions were enhanced to NNLL+NLO by multiplying by an 
appropriate k-factor, and were then combined with ttH 
cross sections from Ref. [32]. For the mass values of the 
heavy scalar considered in this analysis (between 260 and 
320 GeV) the combined cross section from tH and ttH 
reached a value as high as 25 fb at ^fs = 8 TeV. 

The statistical combination was done using the tech¬ 
niques described in Section 3. Firstly, the branching ra¬ 
tios of H —> hh and H —> VV were fixed by minimising a 
y 2 determined from experimental results. These branch¬ 
ing ratios were used as inputs for a combined y 2 , which 
was calculated while floating the free parameters (3 g and 
niH- For each mass point, f3 g was marginalised such that 
the combined y 2 was minimised. Errors on marginalised 
parameters were calculated from identifying the points in 
parameter space which differ by one unit of y 2 above and 
below the minimised value. 
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5. Results and discussion 


In this analysis the global y 2 is minimized for differ¬ 
ent to# hypotheses. The technical part of the analysis 
was done using a scan of mass points, starting at to# = 
260 GeV and going up in 5 GeV steps until 320 GeV. 
Points in between these were reached by an interpolation. 
The other parameters of the model are fixed by a number 
of constraints. Firstly, the branching ratio of H —> hh is 
set to a value that is best fit against the current di-Higgs 
boson resonance search limits set by ATLAS and CMS. 
Secondly, the branching ratio of H —► VV is determined in 
the same way using ATLAS and CMS limits from searches 
for H — > VV at high masses. The remainder of the decay 
of the heavy scalar is assumed to be H —► hyy. Finally, 
the parameter /3 g is constrained by fitting the ATLAS and 
CMS Higgs boson pr spectra, as well as excesses in top as¬ 
sociated Higgs boson production. There may exist other 
decay modes, such as H —► yy. The invisible decay is not 
considered here. Adding in other decay modes would not 
change the final results of the analysis, although it would 
allow us to further constrain the parameter f3 g . 

Calculating and minimising the y 2 described in Sec¬ 
tions 3 and 4, it was found that the lower values of ?n# fit 
the experimental data better than the higher values. Out 
of the mass points considered, the to# = 270 GeV point 
was able to minimise the y 2 value the most. This point 
was determined using the best fit values of the branching 
ratios as BR(7J —> hh) = 0.030 ± 0.037, BR (H —> ZZ) = 
0.025 ± 0.018 and BR {H -a WW) = 0.057 ± 0.041. The 
parameter was best fit at the value of 1.5 ±0.6. The er¬ 
rors on these quantities correspond to a ler deviation from 
the mean value. An indication of this parameter fitting the 
ATLAS and CMS pt spectra can be seen in Fig. 1. The 
fits to the pt spectra were also able to constrain the mass 
of the DM candidate; for to# = 270 GeV, m x was best 
fit at 60 GeV. This is very close to rrih/ 2, which naturally 
leads to the suppression of the branching ratio of h —» yy, 
and it is consistent with current direct search limits. The 
distribution with the filled band corresponds to the pre¬ 
diction made by the SM. The width of the band indicates 
the size of the uncertainties on the ggF process, according 
to the conservative scheme discussed in Section 4. These 
uncertainties are incorporated in the y 2 . The dotted line 
shows the contribution from H —> hyy as well as H ^ hh. 
The solid line corresponds to the sum of the SM and BSM 
components. 

When interpolating between mass points, the combined 
minimised y 2 is found to be smallest at the value to# = 
272 GeV, with upper and lower errors being 12 GeV and 
9 GeV, respectively. This can be seen in Fig. 2 where the 
solid blue line shows the lowest value of the minimised y 2 , 
and the dotted blue lines show a lcr deviation from the 
value. The minimised value of y 2 has a lowest value of 
0.72 per degree of freedom in the fit. It should be noted 
here that when comparing the BSM hypothesis to a null 
hypothesis (the SM with a 125 GeV Higgs boson), the im- 



m H [GeV] 


Figure 2: A scan of minimised \ 2 values as a function of the free 
parameter m#. This was constructed by minimising an additive y 2 
with contributions from all of the experimental results in Table 1. 
These results were compared to the BSM prediction described mostly 
in Section 4. 

provement on explaining experimental data just surpasses 
a 3a effect at the best fit point, as can be seen in Fig. 3. In 
this figure, the large significance around to# = 260 GeV 
can be attributed to the large pp -A H —► hh cross sections 
in most of the ATLAS and CMS di-Higgs boson resonance 
search results. It is also relevant to note that results re¬ 
ported here do not change significantly if the NNLO cor¬ 
rections on ggF±lj discussed in the previous section are 
not applied. 

The most immediate consequence of the phenomeno¬ 
logical model considered here is the appearance of inter¬ 
mediate missing transverse energy in association with h. 
In addition, data appears to display more jets in associa¬ 
tion to h than expected in the SM. This applies both to 
the inclusive production and the production in association 
with top quarks. Enhanced QCD radiation in the produc¬ 
tion of H compared to that of direct h production may 
not be sufficient to explain the effect. As a result, one can 
also consider the decay of the hypothetical intermediate 
particle (discussed in Section 1) into hadronic jets. These 
effects would lead to an enhanced jet multiplicity in the 
inclusive production of h as well as in association with top 
quarks. This also opens the opportunity for a resonant 
structure in the hjj spectrum. 

In the light of this discussion, presented below is a list 
of possible smoking guns for investigation with Run 2 data: 

• Higgs boson in association with moderate and large 
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Figure 3: A scan over mu of the test statistic Xsm — Absm employed 
for quantifying the significance of the effect of the BSM model. The 
la, 2cr and 3cr bands for one degree of freedom are shown in red, 
green and blue, respectively. The SM hypothesis corresponds to the 
existence of the SM Higgs boson, but the complete absence of the 
heavy scalar (i.e. the pp —> H cross section is set to 0 ). 

missing energy. 

• Higgs boson in association with at least two hadronic 
jets, including a resonant structure in the hjj spec¬ 
trum. 

• Higgs boson in association with top quarks and large 
missing energy (greater than 100 GeV) 

• Higgs boson in association with top quarks and at 
least two additional hadronic jets, including a reso¬ 
nant structure in the hjj spectrum. 

• Resonant structure in the VV invariant mass spec¬ 
trum. 

• Resonant structure in the hh invariant mass spec¬ 
trum. 

• Missing energy recoiling against a high transverse 
momentum jet. 

It is beyond the scope of this paper to quantify the 
amount of integrated luminosity required for the ATLAS 
and CMS experiments to declare a 5er effect. However, 
it is noted that the design of dedicated data analysis will 
greatly enhance the sensitivity reached with the available 
results from Run 1 studied here. 

6. Conclusions 

ATLAS and CMS data results comprising measure¬ 
ment of the differential Higgs boson transverse momen¬ 
tum, searches for a di-Higgs boson resonance, the Higgs 


boson in association with top quarks, and VV resonances, 
have been considered against the hypothesis of a heavy 
scalar. The analysis yields a best fit result with a 3er ef¬ 
fect. The hypothetical heavy boson mass is measured to 
be 272±l 2 GeV. In this setup the heavy scalar would pre¬ 
dominantly decay into at least one Higgs boson in associ¬ 
ation with missing energy or hadronic jets. These results 
are obtained on the basis of analyses that are not opti¬ 
mized for the search of the heavy scalar discussed here. It 
is expected that dedicated optimizations will significantly 
enhance the sensitivity to the search. A number of final 
states have been identified that can serve as smoking guns 
for the discovery of the heavy scalar discussed here. 
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